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INTRODUCTION 


Impact  loadings  from  explosive  sources  generally  deform  structural 
components  at  high  strain  rates.  These  strain  rates  are  manifested 
via  elasto-viscoplast ic  structural  response  in  many  materials  such  as 
steel,  titanium,  and  reinforced  concrete.  It  has  been  shown  experi¬ 
mentally  that  the  dynamic  load  carrying  capacities  of  structures  are 
significantly  greater  than  the  corresponding  static  values.  Survey 
papers  of  experimental  work  undertaken  to  study  structural  response 
to  impulsive  loading  have  been  written  by  Jones,  et.  al .  | 27]  and 
Rawlings  [28], 

The  endeavor  to  understand  and  mathematically  model  the  nonlinear, 
inelastic  deformation  process  has  led  to  various  paths  of  development. 
From  an  engineering  perspective,  concern  focuses  on  establishing 
techniques  which  capture  the  essence  of  the  complex  material  behavior 
and  are  suitable  for  the  analysis  of  realistically  modelled  structures. 
Therefore,  this  paper  discusses  the  development  of  an  efficient, 
practical  and  theoretically  sound  method  of  analyzing  the  dynamic 
response  of  structures. 

A  generalization  of  the  viscoplastic  constitutive  model  as 
defined  by  Perzyna  [1]  has  been  employed  in  the  formation  of  a 
system  of  equations  defining  a  viscoplastic  shell  model  in  the 
general  manner  proposed  by  Bieniek  and  runaro  [2]  for  elasto- 
plastic  shells.  This  resulting  shell  model  is  then  incorporated 
into  the  nonlinear  finite  di f fe rence/ f in ite  element  code  EPSA 
(Elasto-Plastic  Shell  Analysis),  ]3],  |4].  The  shell  model 
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developed  in  the  course  of  this  work  utilizes  the  shell  membrane 
strains  and  curvatures  as  the  kinematic  variables  and  the  shell 
stress  resultants  (membrane  forces  and  moments)  as  the  dynamic 
variables.  When  compared  to  the  alternative  through-the-thickness 
integration  method  of  calculating  stress  resultants,  this  approach 
offers  considerable  advantages  in  terms  of  both  computer  storage 
and  speed.  Additionally,  the  model  can  be  fit  to  experimental 
data  representing  various  loading  combinations  (biaxial  bending 
and  stretching)  at  various  loading  rates  on  plate  and  shell 
specimens. 

The  shell  model  approach  can  also  be  applied  to  non-metal  lie 
materials.  In  particular,  a  reinforced  concrete  shell  model  can 
be  formulated,  thereby  avoiding  the  difficulties  involved  in  the 
through-the-thickness  approach  arising  from  stress  re-distribution 
due  to  cracking  and  debonding. 

This  paper  presents  the  qualitative  and  quantitative  aspects 
of  the  viscoplastic  shell  model  through  an  examination  of  a  number 
of  multi-dimensional  loading  combinations  at  various  loading  rates 
for  a  particular  material  type.  The  incorporation  of  the  model 
into  the  EPSA  code  is  outlined,  and  a  study  is  presented  of  the 
effects  of  viscoplastic  material  response  on  the  dynamic  behavior 
of  a  cylindrical  panel  in  vacuo  and  a  closed  cylinder  in  an  acoustic 


medium. 


II  SURVEY  OF  EXISTING  VISCOPLASTIC  THEORIES  AND  THEIR  USE 
IN  SHELL  ANALYSIS 

Continuing  research  into  the  viscoplastic  behavior  of  crystalline 
solids  has  proceeded  from  both  the  microscopic  and  macroscopic 
point  of  view. 

Bingham  [6]  formulated  for  the  case  of  one  dimensional  shear 
a  relationship  represented  by  Figure(l).  Thic  concept  was  extended 
to  three  dimensional  states  of  stress  by  Hohenemser  and  Prager  [7] 
with  the  resulting  expression  for  the  plastic  strain  rate  £  P 

ij 


in  the  form 


•  P  _  _L  *  2  __ 

^  2y  ^rrr 


where  is  the  second  invariant  of  the  deviator ic  stress  4  . 

(U  and  k  are  material  constants) . 

A  generalization  of  these  basic  ideas  has  been  developed  by 
Perzyna  [1].  Included  in  the  specification  is  the  existence  of  a 
quasi-static  yield  function 

f(0  ,  ej  ) 

F  =  - 3J - H  _  1  (2) 

K 

p 

where  £ .  . 

lj  P 

K  *  <(Wp)  -  k  atj  de..  (3) 

o 

is  the  work  hardening  parameter  which  describes  how  the  quasi¬ 
static  yield  surface  deforms  during  the  inelastic  process. 

The  dependence  of  the  plastic  strain  rate  tensor  on 
the  stress  intensity  is  given  by 

3^  (1> 


where  y  and  [ F ]  are  material  response  functions.  The  direction  of 


the  viscoplastic  strain  rate  vector  is  directed  along  the  normal 
to  the  subsequent  dynamic  yield  surface. 

The  dependence  of  the  yield  criteria  on  the  strain  rate  and 
the  dynamic  yield  criteria  for  elasto-viscoplastic  work-hardening 


materials  is  given  as 

f(a..,  e.?)  =  <(w  ) 

iJ  P 

where 


r  n  1/2 

(I  p) 

v  2  ;  JL  9f  3f 

Y  2  3V  9V 


1 

2 


P 

ij 


(5) 


(6) 


Figure (2)  shows  a  dynamic  stress-strain  relation  for  elastic 
viscoplastic  work-hardening  material  in  simple  tension. 

A  mathematical  description  of  the  inelastic  metal  deformation 
process  proposed  by  Cernocky  and  Krempl  [8]  is  based  upon  the 


following  relation  for  the  plastic  strain  rates 

•  P  =  r  gkJl  Gk£ 

ij  ^ijkS,  K[V] 

where  is  the  elastic  compliance  tensor,  is  the  stress 

tensor  corresponding  to  the  "equilibrium  stress-strain  curve",  T 


is  the  second  invariant  of  (o^  -  and  <[T]  is  a  material 

function. 


Motivated  by  dislocation  dynamics,  Bodner  and  Partom[9] 
proposed  a  functional  relationship  between  the  second  invariant 
of  the  plastic  strain  rate  and  the  second  invariant  of  the 
deviatoric  stress 

°2  =  FfJ2]  (8) 

where  the  form  of  the  functional  is  chosen  to  reproduce  the 
stress-strain  curves  of  a  particular  material. 

Werne  and  Kelly  [ 10]  arrived  at  a  model  of  metal 
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viscoplasticity  from  dislocation  dynamics  considerations 

•  8bo  (9> 

where  $  is  a  nonlinear  function  of  the  deviatoric  stress  invariant 

H  and  £  are  defined  in  terms  of  their  evolutionary  equations  and 

represent  the  dislocation  density  and  the  mobile  fraction, 

respectively,  is  Burger's  vector.  This  model  was  shown  capable  o 

reproducing  various  types  of  uniaxial  stress-strain  curves. 

Material  models  incorporating  viscoelastic  behavior  have  also 
been  introduced.  Naghdi  and  Murch  [25]  have  postulated  a  theory  of 
viscoelastic/plastic  solids  which  reduces  to  that  of  linear 
viscoelasticity  and  to  that  of  classical  plasticity  in  limiting 
cases. 

Inelastic  analysis  includes,  in  addition  to  a  mathematical 
description  of  the  metal  constitutive  behavior,  a  solution 
scheme  for  the  governing  nonlinear  equations  of  motion. 

Zirin  and  Krempl  [11]  used  a  finite  element  technique 
to  solve  plane  stress,  plane  strain  and  axisymmetrical  problems 
incorporating  Krempl's  [  8]  viscoplastic  model.  A  particular 
problem  investigated  was  a  thick-walled  axisymmetrical  cylinder 
under  internal  pressure. 

Werne  and  Kelly  [10]  used  the  HONDO  code  [12] 
response  of  a  cylinder  and  a  bar  in  uniaxial  tension. 

A  numerical  analysis  of  the  elastic  viscoplastic  response 
of  an  axisymmetrical  spherical  shell  was  performed  by  Takezono, 
et.  al.  [13]  using  Perzyna's  constitutive  model. 

The  incorporation  of  the  viscoplastic  shell  model  into 


the  EPSA  code  for  usage  in  solving  three  dimensional  shell  type 
structural  problems  required  the  selection  of  a  general  form  for 
the  viscoplastic  constitutive  relations.  The  representation  chosen 
was  based  upon  Perzyna's  model  [1],  The  governing  arguments  for 
this  selection  were: 

1.  The  functional  form  of  the  constitutive  relations  is  of 
sufficient  generality  so  that  a  broad  range  of  material 
behavioral  patterns  could  be  represented. 

2.  Experimental  verification  by  Hayashi  and  Tanimoto[26] 

of  a  Perzyna  type  formulation  for  describing  the  response 
of  annealed  aluminum  to  impulsive  torque  and  tension 
loadings . 

3.  Retention  by  the  constitutive  model  of  certain  basic 
elements  of  the  classical  theory  of  plasticity  which  have 
proved  to  be  a  practical  and  accurate  tool  of  analysis 
for  elasto-plastic  structures.  In  particular,  uniqueness 
arguments  of  the  type  made  by  Drucker[24]  are  applicable. 

4.  The  model  incorporates  multi-dimensional  stress-strain 
states  including  loading  and  unloading  paths. 

5.  A  shell  model  formulation  could  be  established 
incorporating  the  material  constitutive  equations. 

6.  A  stable  solution  scheme  for  the  viscoplastic  shell 
equations  could  be  established. 

7.  The  model  is  effective  in  capturing  complex  material 
behavior  while  still  being  cost  effective  from  a 
computational  point  of  view. 


MATERIAL  CONSTITUTIVE  EQUATIONS 


Perzyna  postulated  a  set  of  constitutive  equations  to 
represent  viscoplastic  strain  hardening  behavior  for  arbitrary 
loading  histories  [[1].  Included  in  this  specification  is  the 
direction  of  the  viscoplastic  strain  rate  and  the  magnitude  of 
allowable  stress  beyond  the  static  yield  surface. 

The  particular  case  of  constitutive  relations  chosen  is 
that  of  elastic  visco-perf ectly  plastic  material  behavior 
for  which  the  static  yield  function  in  stress  space  is  the  von 
Mises  yield  condition 

F(.  >  -  1  *i-1  4i-1  _  , 

F(<5ij}  2  r2  1 

1  *  o 

where  4,,  are  the  deviatoric  stresses,  K  =  -=—  and  a  is  the 

±3  yp 

uniaxial  static  yield  stress. 


The  rate  of  increase  of  the  inelastic  components  of  the 
strain  tensor  is  a  function  of  the  excess  stresses  beyond  the 
static  yield  surface (overstress) .  The  direction  of  the  viscoplastic 
strain  rate  vector  is  along  the  normal  to  the  subsequent  loading 


surface.  These  relations  are  mathematically  stated  as 

-  YK»[F(i  )1  —  -A  —  (11) 

3  3  Wlj  ^IJ 

where  y  and  $[F]  are  material  functions  which  can  be  chosen  to 
represent  the  results  of  experimental  tests  on  the  dynamic  behavior 
of  a  particular  material.  Current  available  experimental  methods 
consist  of  one-dimensional  stress-strain  tests  performed  using  a 
hydropneumatic  machine  or  a  hopkinson  bar  technique  [14]. 

Figure (3)  shows  a  one  dimensional  yield  stress  versus  strain 
rate  curve.  Manjoine[17]  states  that  this  "S"  shaped  curve  is 
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1 


typical  of  metals  tested  within  a  given  range  of  strain 

rates.  The  particulars  of  the  curve  vary  for  each  material  type. 

An  approximation  or  fit  to  this  characteristic  yield 

stress  versus  strain  rate  curve  is  made  by  the  following  expression 

n2 

P  I/O.  a  (F) Z 

tll  “  al(F)  +  n„  (12) 


where 


ft  =  (-Z-)  -  1 

l  a 

o 


and  a.  is  a  limiting  one  dimensional  yield  stress  achieved  from 
dynamic  tests  at  the  highest  of  strain  rates  of  interest. 

Note  that  equation  (12)  can  be  reduced  to  either  a  convex 
or  concave  power  law  (with  or  without  an  upper  limiting  stress) 
representation  of  the  yield  stress  versus  strain  rate  behavior. 

The  power  law  representation  has  been  prevalent  in  previous  works 
[15]  and  [16]. 

This  one  dimensional  experimental  data  can  be  represented  in 

Perzyna's  formulation  by  assuming  9  [F]  =  F  and  using  a  suitable 

functional  form  for  Y  .  Then,  specializing  equation  (11)  to 

2 

one  dimensional  behavior  a  a 

.  P  2  11  ,11  .... 

en  vs  y~  (~i  -  “  U4) 

0  °o 

or  in  terms  of  the  static  yield  function  F(^  ) 

hi  vlYF(VF  + 1}  (i5) 

Equations  (12)  and  (15)  state  that  the  material  response 


function  Y  is 


VST _ 

2  I  1-1/n 


'F  +  1 


a2  (F) 


(fl  -  F)  VF  +  1 


K 
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Appropriate  values  of  the  material  constants  a^  ,  a^  , 
nl  *  n2  an£*  n3  are  ^eterm^ne^  by  curve  fitting  to  one  dimensional 
stress-strain  rate  experimental  data. 

As  an  illustrative  example,  the  uniaxial  behavior  of  a  mild 

0  3 

steel  within  the  strain  rate  range  of  (10)  to  (10)  1/sec  is 
investigated. 

The  experimental  data  of  Clark  and  Duwez  [18]  is  fit  to  via 
the  parameter  values  shown  in  Figure (4). 

The  general  character  of  equation  (16)  is  that  Y  approaches 
infinity  at  the  asymptotes  F  =  0  and  F  =  F^  .  This  behavior  can 
also  be  represented  by  the  expression 

ai _ 

Y  *  — n  n~  (17) 

1  2 

F  X(Fl  -  F)  1 

The  curve  fit  and  parameter  values  using  equation(17) 


are  shown  in  Figure  (5)  for  the  mild  steel  material. 


IV 


SHELL  CONSTITUTIVE  EQUATIONS 


Shell  structures  resist  external  loading  by  developing 
internal  biaxial  forces  and  moments.  Two  alternative  computational 
approaches  for  obtaining  these  internal  stress  resultants  are 
presented. 

a.  Through-The-Thickness  Integration 

For  this  technique,  the  time  increments  of  strains  at  any 
point  of  the  shell  are  expressed  in  terms  of  the  strain  and 
curvature  increments  of  the  middle  shell  surface. 


Vu<2) 


Vu  +  v-  • z 


'il  ij 

The  time  increments  of  shell  forces  and  moments (stress 


(18) 


resultants)  are  determined  by 
h/2 


h/2 


*nn 


ij 


A^O.,  ^  dz 


-h/2 


'ij 


Vij 


V 


-h/2 


ij 


z  dz 


(19) 


The  integrals  are  numerically  computed  by  dividing  the  shell 
thickness  into  "k"  layers  and  assuming  a  linear  stress  distribution 
within  each  layer. 

The  stress  increments  are  computed  from  the  material 

constitutive  equations 


V 


ij  *  °ij  "  °ij  *C  ijk*  AN(£k<t  ‘  £kJ^ 


(20) 


where  O  is  the  stress  tensor  after  "N"  time  increments, 

C.,,  ,  is  the  elastic  moduli  matrix.  The  total  strain  increment 
ijkl 

to  be  separated  into  elastic  and  inelastic  components. 


Equation  (20)  may  be  interpreted  as 


0  N  *  O  E1Tr  -  c  A  £  p  (22) 

ij  ij  ljkJi  TTkA 

ElTr 

where  is  the  stress  tensor  based  upon  a  trial  assumption 

th 

of  total  elastic  behavior  occuring  within  the  N  time  step. 

°ijE1Tr  “  °ij  1  +  Cijk£  Vkj£  (23) 

Writing  the  viscoplastic  relations  (equation(ll) )  in  time 


incremental  fashion 

- 


Vii  ■  4tY  KF^ii  > 


N  3F  9F 

ij  '  3*  N 

ij  ij 


X  is  a  non-negative  variable  which  must  be  determined  as  follows; 

A  -  -  C  «5) 


Noting  N  =  2  and  *  2G  4  Equation  (25)  becomes 

3*ij  K  J 

A  -  YAtKT(41JE1Tr  -  ^  Ajj  )  (26) 

K 

Equating  arguments  of  function  FCa^j)  yields, 


1+^A 

K 


X  =  yAtK  (- 


1.  ElTr  ElTr 
2  *il  *11 

-7G  ^  2 

(1  +  ^fX)  k 


-  1)  (28) 


Rearranging  equation  (28)  gives 

C^J)  X3  +  (  J  x2  +  (1  +  ^S^t)  X 

K  K  K  K 

(29) 

“YAtK  (F(4ljE1Tr))  =  0 

The  real  root  is  extracted  from  a  cubic  equation  solver  of 


X  [  v  ,  A  t ,  K  ,  G  , F (a 


)  ]  =  0  where  Y  is  given  as  described 


in  Section  III. 
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vT1-"  V:  Sr  A  4 . . ■  a  AA. ■  • V  *  v  -  A  •  A  .  • 


■, 1  ^  '  u  i 


1  •  ■ 


and 


»0  ■  V 


M. 


V  /6 


where  a  is  the  static  yield  strength, 
o 


The  quantities  M^  ,  which  will  be  referred  to  as 


!l  hardening  parameters  "  are  defined  by 


M 


If  F  =  1  and  N.  .  +  M.  .  >0:  dM.*  =  3.  .  (1-FT  )  —  3k.  ,P 

3Nij  l 2  9Mij  ij  ij  ij  L  Kq  ij 


If  F  <  0  or  N.  .  +  -r— •  M.  .  <0:  dM.*  =  0 

3N. .  ij  dM. .  ij  —  ij 


(36) 


ij 


ij 


The  elastic  law 


Sij  =  EijK£(ekS,  “  eU) 


(37) 


is  used  where  E. ,  is  the  elastic  shell  stiffness  matrix, 
ijkl 


The  above  quasistatic  shell  formulation  is  derived  from  that 
proposed  by  Bieniek  and  Funaro  [2]  for  elasto-plast ic  shell  analysis 
Alterations  have  been  made  in  the  specification  of  the  variable  yield 
surface,  the  limit  surface,  and  the  hardening  law,  resulting  in  a  more 
accurate  representation  of  the  quasistatic  elasto-plast ic  shell  behavior. 
In  order  to  generalize  the  model  to  include  strain  rate  effects  along 
the  lines  of  equation  (11),  the  plastic  strain  rate  tensor  is  assumed 
to  depend  on  the  stress  resultant  intensity  through  the  associated  flow 
rule 


ei.P  =  YRK<J>[F(Si.)] 


3F(S  . )  3F (S  .  . ) 

- !L_  =  X - JUL_ 

as. . 

i.i  ij 


(38) 


Recalling  the  work  of  section  III,  one  obtains 


ElTr 


3F 


A  =  VtKF(sij  -Eiikf  A;^:) 


(39) 


kit 
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Because  equation  (39)  cannot  be  solved  in  closed  form  for  A 


the  equivalent  form 

MYr,  At,  K,  Eljfcr  F(SiJE1Tr)>3|^  ]  -  0  (40) 

is  solved  iteratively  using  a  modified  regula  falsi  method, 
is  a  physical  function  which  should  be  chosen  to 
represent  the  results  of  multiaxial  tests  on  plate  and  shell 
specimens  of  a  particular  material.  These  experiments 
should  consist  of  biaxial  bending  and  stretching  tests  conducted 
at  various  loading  rates  of  concern.  The  present  lack  of  such 
experimental  data  precludes  the  possibility  of  choosing  this 
function  to  fit  measured  results.  Therefore,  at  this  time,  Y 

K 

will  be  chosen  such  that  the  produced  stress  resultants  obtained 
via  the  proposed  viscoplastic  shell  model  satisfactorily  represent 
those  obtained  using  the  through-the-thickness  approach  for 
various  loading  paths. 

For  the  mild  steel  material  [18],  previously  discussed,  the 

proposed  form  of  the  function  Y„  is 

R 


VV  XM  1 


where 


and 


ni  an  2  n7 

(Fr)  (1  -  (32)  Fr) 


N  “2“3 


*  I„  +  A„A 

A1^IN»  J  *  Al^'  *  ' 

h  +  Vm 

C1(IN  +  XN  XN  *  +  C2(IM  )  +  °3  I XNM 


(41) 

(42) 

(43) 


This  formulation  respresents  an  eight  parameter  model 

n^  ,  n^  ,  A^  ,  A^  ,  A^  ,  )  whose  values  are 

chosen  such  that  the  viscoplastic  shell  model  results  correspond  to 


the  calculated  through-the-thickness  viscoplastic  shell  behavior 
for  various  loading  rates  and  paths  for  the  mild  steel  specimens. 
A  discussion  of  the  determination  of  these  parameters  will  be 
given  in  Section  V. 


COMPARISON  OF  THROUGH-THE-THICKNESS  AND  SHELL  MODEL 


RESULTS  FOR  THE  BEHAVIOR  OF  MILD  STEEL 


Multi-dimensional  structural  response  involves  the  determination 

of  biaxial  shell  forces  and  moments  for  various  combinations  of 

induced  axial  and  bending  strains.  With  this  is  mind,  a  series 

of  loading  cases  is  introduced  ,  with  each  case  reflecting  a 

particular  stress  state  feature.  Figure  (6)'  presents  a  listing 

of  particular  loading  cases  considered. 

Continuing  the  illustrative  example  of  mild  steel  behavior, 

the  viscoplastic  shell  model  parameters  of  equations  (42)  and  (43) 

were  determined  such  that  the  shell  model  results  matched  the 

through-the-thickness  results  for  the  loading  paths  listed  in 
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Figure  (6)  for  loading  rates  (10)  to  (10)  1/sec. 

For  the  mild  steel  considered,  the  determined  material 
constants  are 


«  3/4 

n2  *  1/4 

20.0 

^2  *  1.0 

a3  *  1.0 

1.0 

C2  -  0.4 

C3  =  0.8 

For  the  loading  cases  analyzed,  various  features  of  the 
viscoplastic  shell  behavior  are  presented  in  Figures (7)  to  (23). 

Figure  (7)  shows  the  variability  in  the  axial  force  as  a 
function  of  the  applied  strain  rate.  An  increase  of  the  yield 
force  is  seen  for  increasing  strain  rates  and  a  corresponding 
reduction  in  yield  force  for  decreasing  strain  rates.  It  is  seen 
that  a  rise  in  the  strain  rate  raises  the  flow  stress  rapidly 
from  its  previous  stationary  value  and  as  straining  continues 
at  the  more  rapid  rate,  the  stress  approaches  the  value  it 


■t 


would  have  had  if  the  entire  straining  process  was  performed  at 
this  rapid  strain  rate.  Thus,  the  effects  of  strain  rate  history 
are  erased.  Experimentally,  it  has  been  shown  [19]  that  for 
mild  steel,  strain  rate  history  effects  are  minimal  at  ambient 
temperatures . 

For  the  viscoplastic  shell  behavior  in  one  dimensional  bending. 
Figure  (8)  presents  the  moment-curvature  diagrams  for  various 
loading  rates.  Figure  (9)  shows  the  increase  of  the  limit  force 
and  the  limit  moment  due  to  strain  rate  effects.  One  notices  in 
Figure  (9)  that,  at  a  particular  loading  rate,  the  limit  moment 
is  not  as  dynamically  enhanced  as  the  limit  force.  This  is  due 
to  the  variability  in  the  strain  rates  through-the-thickness  of 
the  shell  when  deforming  in  a  bending  mode. 

Loading  cases  three  through  six  represent  a  variety  of 

biaxial  stress  resultant  states  dominated  by  a  particular  force 

and/or  moment.  Various  radial  loading  paths  of  stretching  (Ae)  and 

bending  (Ak)  strain  increments  were  prescribed  to  produce  the 

biaxial  stress  resultant  states  =  N  ,  ^  =  aN  ,  =  M  and 

M2  =  3m  .  These  loading  cases  were  investigated  for  loading 
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rates  varying  from  (10)  to  (10)  1/sec.  For  the  sake  of  brevity, 
the  loading  rate  of  120  1/sec  is  chosen  to  illustrate  the 
shell  model's  applicability.  Each  of  the  featured  cases  consists 
of  a  cyclic  loading  curve.  The  viscoplastic  shell  model  captures  the 
essential  features  and  characteristics  of  the  through-the-thickness 


integration  calculations,  as  shown  in  Figures  (10)  to  (14)  . 
Figures  (15)  to  (18)  illustrate  the  variability  in  the 


stress  resultant  limit  state  for  the  condition  of  equal  biaxial 
forces  and  equal  biaxial  moments.  The  ratio  of  stretching  (Ae)  to 
bending  (Ak)  is  varied  to  produce  various  combinations  of 
forces  and  moments.  It  is  clear  that  the  viscoplastic  shell 
model  produces  a  very  good  definition  of  the  effect  of  strain  rates 
on  the  dominating  stress  resultant  quantities.  The  existing  small 
deviations  are  seen  to  be  independent  of  the  loading  rate  and  are 
therefore  due  to  approximations  inherent  in  the  specification  of  the 
quasi-static  yield  surface. 

Loading  cases  eight  through  twelve  represent  equal  absolute 
magnitude  biaxial  forces  and  biaxial  moments,  with  differing  signs. 
Figures  (19)  through  (23)  illustrate  these  cases  and  show  the  good 
agreement  which  exists. 

Therefore,  for  the  various  loading  combinations  which  can 
occur  in  shell  structures,  it  is  shown  that  the 
model  provides  a  very  good  definition  of  the  viscoplastic  shell 
behavior  as  determined  by  the  through-the-thickness  integration 
procedure.  The  selection  of  suitable  but  different  model  parameters 
should  provide  the  same  favorable  approximations  for  other 
structural  metals. 

When  biaxial  experimental  data  becomes  available,  the 
model  parameters  can  be  fit  to  the  data.  It  has  been  shown  that 
for  dynamic  one  dimensional  bending,  deviations  do  exist  between 
measured  moments  and  those  computed  using  the  through-the- 
thickness  integration  technique [ 20 ] . 


VI 


INCORPORATION  OF  VISCOPLATIC  MODEL  INTO  EPSA  CODE 


The  EPSA  [3  ]  code  has  been  developed  for  the  dynamic 
response  analysis  of  stuctures  in  an  acoustic  medium  , 
including  both  plastic  collapse  and/or  local  instability 
dynamic  buckling)  of  the  structure.  The  following  is  a  brief 
description  of  the  theoretical  basis  of  the  code: 

The  structural  equations  of  motion  are  derived  for  the 


principle  of  virtual  work 


' 


{S}  {6e}dR  - 


{P}  {6u}  dR  + 


jp{u}{6u 


}  dR  =  0  (4A) 


With  (u)  =  (Ul,  u2,  w)T 

A  finite  difference  /  finite  element  method  is  used  to 

discretize  the  non-linear  equations  of  motion.  The  surface  of 
the  region  R  is  covered  with  a  quadrilateral  mesh  of  "j"  elements, 
each  element  of  area  A^. 

Each  arbitrarily  shaped  quadrilateral  shell  element  is 
defined  by  four  corner  nodes,  with  each  node  having  three 
translational  and  no  rotational  degrees  of  freedom.  Spacial 
derivatives  are  expressed  in  terms  of  discrete  nodal  displacements 
via  an  irregular  finite  difference  technique.  A  two  dimensional 
Taylor  series  expansion  in  irregularly  shaped  meshes  is 
employed  [5]i. 

The  integrals  over  region  R  of  equation  (44)  are  replaced  by 
sums  of  integrals  over  A^  to  obtain  the  following  system  of 
ordinary  differential  equations: 


j 

[M]  {q}  -  l  [{PL  ~  {F}  ] 
i-1 


(*5) 
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where  {q}  is  the  nodal  displacement  vector  for  the  structure, 

[M]  is  the  lumped  mass  matrix  ,{P}  represents  the  vector  of 
external  forces  acting  on  the  nodes  of  the  structure,  {F}  is 
the  vector  of  equivalent  internal  grid  point  forces. 

The  system  of  equations  for  the  nodal  displacements  is 
integrated  in  time  using  an  explicit  central  difference  scheme. 

The  Donnell-Vlasov  nonlinear  kinematic  equations  of  shell 
theory  are  used.  The  geometric  nonlinearities  are  accounted  for 
in  the  strain  -  displacement  relations 

{e}  *  [B]  {q}  (46) 

The  viscoplastic  shell  equations  relate  tne  stress  resultant 
vector  to  the  shell  strain  rate  vector.  These  equations  are 
solved  at  each  time  step.  The  magnitude  of  the  time  step  is 
governed  by  the  Courant  criterion  for  the  integration  of 
the  spatially  discretized  equations  of  motion. 

The  use  of  the  viscoplastic  shell  theory  results  in  an 
increase  in  the  efficiency  of  the  computational  procedure  as 
compared  to  the  through-the-thickness  integration  technique. 

This  occurs  because  the  shell  theory  employs  nine  memory 
parameters  per  element  versus  the  42  (  6  *  Number  of  Layers) 
quantities  for  the  through-the-thickness  technique.  Also,  the 
solution  procedure  for  the  shell  model  requires  approximately 
50%  less  computer  time  than  the  through-the-thickness  technique. 


VII  VISCOPLASTIC  SHELL  ANALYSIS  USING  EPSA 


REPRESENTATIVE  EXAMPLES 

The  EPSA  code  with  the  inclusion  of  a  viscoplastic  shell  model 
was  used  to  investigate  the  elastic  viscoplastic,  large  deflection, 
transient  response  of  shell  structures. 

An  analysis  was  performed  to  determine  the  response  of  a 
clamped-edge,  mild  steel,  cylindrical  panel  subjected  to  impulsive 
loading  by  the  sheet  explosive  loading  technique.  The  geometric 
properties  of  the  panel  are  shown  in  Figure  (24). 

Analysis  of  the  structure  by  EPSA  consisted  of  a  one 
hundred  element  mesh  configuration  for  the  quarter  model, 

Figure  (25).  Initial  conditions  consist  of  initial  radial  velocities 
obtained  by  equating  the  impulse  imparted  by  the  detonated  sheet 
to  the  total  impulse  experienced  by  the  structure. 

Figures  (26)  and  (27)  illustrate  the  results  of  two  EPSA 
analyses,  one  using  an  elasto-plastic  shell  model  (quasi-static 
yield  strength),  the  other  using  a  viscoplastic  shell  model.  The 
viscoplastic  material  effects  are  shown  to  reduce  the  permanent 
deflection  at  the  center  of  the  specimen  by  twenty  five  percent. 

The  peak  inner  fiber  circumferential  strain  at  the  mid-point  of 
the  specimen  is  reduced  from  3  percent  to  2  percent  due  to 
viscoplastic  shell  behavior. 

The  second  example  investigated  involved  a  stiffened 
cylindrical  shell  immersed  in  an  acoustic  medium  subjected  to 
shock  loading.  Figure  (28).  The  pressure  function  contains  an 
exponential  decay  in  time  and  a  linear  spatial  decay. 

The  EPSA  quarter  model  of  the  shell  and  internal  stiffeners 
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consists  of  a  1440  element  mesh.  The  fluid/structure  interface  is 
accounted  for  by  means  of  the  Doubly  Asymptotic  Approximation  f 231 . 

Again,  two  distinct  EPSA  analyses  were  performe  to  isolate 
the  effects  of  strain  rate  behavior  upon  the  calculated  response 
of  the  structure.  Figure  (29)  shows  the  calculated  deflection 
pattern  across  the  length  of  the  shell  at  the  circumferential 
zero  degree  line  (  initial  load  impact  position)  fot  the  two 
analyses.  The  variability  in  the  inner  fiber  circumferential 
strains  at  Frame  5  1/2  at  zero  degrees  and  180  degrees  is  shown  in 
Figure  (30)  for  the  two  analyses. 


VIII  CONCLUSIONS 


An  analysis  method  is  presented  for  obtaining  the  dynamic 
response  of  elasto^viscoplastic  structures.  This  technique  consists 
of  the  development  of  a  viscoplastic  shell  theory  and  its 
incorporation  into  the  EPSA  code.  The  effects  of  viscoplastic 
material  behavior  upon  the  response  of  structures  in  vacuo  or 
in  an  acoustic  medium  when  subjected  to  transient  loadings  is 
demonstrated. 

Viscoplastic  material  response  is  shown  to  significantly 
increase  the  resistance  of  shells  to  dynamic  loadings. 

Continuing  development  work  can  encompass  the  following: 

1.  Employ  the  proposed  method  to  analyze  the  viscoplastic 
response  of  mild  steel  frame,  plate  and  shell  specimens 
and  compare  to  available  experimental  data,  [21  ],[22  ]. 

2.  Fit  the  shell  model  parameters  to  additional  material 
types,  i.e.,  titanium  and  high  strength  steels. 

3.  Incorporate  strain  hardening  into  the  material  constitutive 
equations . 

4.  Explore  the  possibility  of  obtaining  experimental 
biaxial  stress  resultants  versus  strain  rate  data  so 
that  the  viscoplastic  shell  model  parameters  can  be  fit 
to  such  data  rather  than  to  the  through-the-thickness 


results . 
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